Transport in Graphene: Ballistic, Diffusive, or Anomalous? 
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We investigate the transport of electrons in disordered and clean graphene devices. We consider a 
geometry where the graphene flake is contacted by narrow metallic leads. Plotting the conductance 
as a function of the position of one of the leads, we can approximate the probability density function 
of the charge flow at the edge which is used to analyze the transport properties with increasing length 
of the device. Moreover, we simulate scanning probe microscopy (SPM) measurements for the same 
devices, which can be seen as a measure for the flow of charge inside the device, thus complementing 
the transport calculations. We compare our analysis to theory describing transport in clean and 
disordered systems. 



I. INTRODUCTION 

As predicted by early band theory studies [1 , charge 
carriers in graphene are chiral quasiparticles that have a 
linear Dirac-like dispersion relation resulting in fascinat- 
ing electronic [2^ and transport properties . As a result 
of the linear energy- momentum relationship, the valence 
and conduction bands intersect at the Dirac points lo- 
cated at the Brillouin zone corners. For neutral graphene, 
the chemical potential lies exactly at the Dirac points. 
Use of a gate potential can shift the chemical potential 
away from the Dirac point and thus tune the charge 
density of the graphene device over a large range (as 
compared to other electronic materials). These tunabil- 
ity of electron density and the linear dispersion are fea- 
tures that distinguish graphene from the two-dimensional 
electron gas in semiconductor heterostructures and other 
metallic conductors. Understanding the nature of charge 
carrier transport in graphene in the low carrier density 
regime (in this linear dispersion regime) and in the pres- 
ence of disorder has received considerable attention [2]-[4] . 

A wide variety of fabrication techniques [5HTQ] has re- 
sulted in graphene devices that operate in different dis- 
order regimes. The mobility of the devices gives a way 
to quantify the level of disorder the electrons encounter. 
The standard method of mechanical exfoliation and de- 
positing graphene flakes on Si02 substrates [5^ can re- 
sult in good quality devices with mobilities approaching 
a few tens-of-thousand cm?V~^ ^ values comparable to 
those found in conventional semiconductor devices. Yet, 
the substrate provides a source of impurities [11] and 
corrugations for the electrons to scatter, limiting the 
mean free path to about a hundred nanometers. Etch- 
ing of the substrate to produce suspended graphene de- 
vices followed by current annealing to remove a substan- 
tial number of impurities increases the carrier mobility 
to values exceeding 200,000 cm?V~^s~^ and mean free 
paths / approaching the size of the device [SI H] . Similar 
sample quality is obtained in devices where the substrate 



is single-crystal hexagonal boron nitride with mobilities 
/i reported as high as 80,000 cm?V~^ jlOj . 

It is apparent that replacing oxide based substrates 
by high-quality crystal substrates or suspending the 
graphene results in devices with improved mobility and 
reduced electron scattering times jSHTOj. In an ideal 
(no impurities or defects) strip of graphene, the mean 
free path is longer than the size of the device result- 
ing in ballistic conductivity. In the case of transport 
at the Dirac point, where the density of states vanishes, 
the conductivity has a minimum value (4e^/7r/i), even 
for ideal graphene [13 . This unusual result stems from 
the propagation of charge carriers through evanescent 
waves [13] . Experimentally, the phase coherent transport 
in graphene devices with short- and- wide geometries has 
been characterized as ballistic [HI [15] indicating evanes- 
cent wave transport. Thus, transport in graphene quali- 
tatively changes from metallic diffusive in Si02 to what 
has been characterized as ballistic in suspended systems. 

The variety of transport mechanisms makes it essen- 
tial to develop measures to clearly distinguish and classify 
the different transport regimes. We begin with a simple 
question - How does the conductance change with re- 
spect to the vertical displacement of the two leads? In 
particular, we study the transport properties of clean and 
disordered graphene devices contacted by narrow metal- 
lic leads where one of the leads can be moved along the 
edge of the device. We propose that this method can 
be used to unveil the mechanisms of electronic transport 
given that the probability density function (PDF) of the 
charge flow at the edge of the device can be approximated 
from the dependence of the conductance to the position 
of one of the leads in the device. Moreover, we simulate 
scanning probe microscopy (SPM) [16 measurements for 
the same devices (but with fixed leads), which sheds light 
on the charge flow inside the device, complementing the 
transport calculations. 
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II. TRANSPORT MEASURES 
A. Shot noise and Fano factor 

Shot noise measurements can determine the statis- 
tics relevant for transport in mesoscopic conductors [17 . 
Shot noise is a consequence of charge quantization and 
can be assessed by the Fano factor defined as the ra- 
tio of noise power and mean current, T = ^n(l — 
^n)/ Z]n=i From the above definition, the Fano fac- 
tor should be zero for conventional ballistic systems, i.e. 
if perfect transmission is present. It has also been estab- 
lished that for Poisson process T = \ and for diffusive 
metallic conductors = 1/3 [18 . 

However, in the case of transport in ideal graphene, cal- 
culations predict that the Fano factor at the Dirac point 
is also 1/3 [13]|19]. A slightly smaller and universal value 
[T — 0.295) was found numerically for disordered sys- 
tems [20l [21] . When the disorder is mainly from rough- 
ness in the substrate, the Fano factor is nearly constant 
but non-universal, above the 1/3 value, and increasing 
with disorder [22 . Transport along an n-p (electron rich/ 
hole rich) junction is selective of those quasiparticles that 
approach the n-p interface almost perpendicularly and re- 
sults in a Fano factor of = 1 — y^1/2 [23 , numerically 
very close to the 1/3 value. 

The Fano factor measured in clean devices is close 
to the value of 1/3 found analytically for ballistic sam- 
ples [l^ and measurements in disordered devices [24 fol- 
low the numerical trends of Ref. [22]. However, these val- 
ues are numerically near each other making the Fano shot 
noise a cumbersome measure of the transport mechanism 
in graphene devices [3]. Therefore, developing an alterna- 
tive quantitative understanding of the carrier dynamics 
is crucial for testing the transport regime of graphene 
devices. 



B. Probability density function 

The transition from diffusive to ballistic transport can 
be quantified by statistical means. Diffusive transport 
in disordered systems is based upon multiple scattering 
of the charge carriers off impurities and can be described 
classically by a random walk, where the mean-square dis- 
placement reads {x{t)'^) ex t", with exponent a = 1. The 
step lengths are described by a Gaussian probability dis- 
tribution (which can be depicted, e.g., using a Galton 
board). Ballistic diffusion corresponds to a = 2, while 
for 1 < a < 2 superdiffusion, also known as a Levy walk, 
is encountered. Levy walks belong to generalized random 
walks in which the step lengths are described by a proba- 
bility distribution with power law tails. In this case, the 
step lengths distribution features an algebraic tail and 
describe stochastic processes that are scale-free. Theoret- 



ical concepts propose Levy walks in solid state systems, 
e.g. in semiconductor antidot lattices [25' and Josephson 
junctions [26]. However, except for amorphous materials 
[271429] , the possibility to observe and study Levy trans- 
port in solid state devices is exceedingly rare. For exam- 
ple, experimental work on heat, sound, and light diffusion 
is generally limited to normal diffusion. Recently, a Levy 
walk for light has been realized by engineering an optical 
material [30] . 

III. COMPUTATIONAL METHOD 

We use the nonequilibrium Green's function (NEGF) 
formalism [31 to calculate the transport properties of a 
normal-conductor / graphene / normal-conductor junction, 
given that in experiments the electronic contacts are usu- 
ally metallic. The metallic lead/graphene interface we 
have investigated follows that of Robinson and Schome- 
rus [32 . An schematic of the geometry of the devices 
investigated is shown in Fig. [l] where rectangular (width 
W and length L) graphene fiakes have armchair bound- 
aries along their longitudinal direction and narrow metal- 
lic leads connected along the zigzag boundary. The two 
narrow leads of width are modeled by semi-infinite 
square lattice regions. 

The tight-binding model Hamiltonian of the device is 
given by 

H = Y^ eic\ci - -fijclcj, (1) 
where is the on-site potential of the lattice site z, q(cJ) 
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FIG. 1: Schematic setup for the transport calculation. 
A rectangular graphene flake, depicted by red and blue 
atomic sites arranged in a honeycomb lattice, is contacted 
between two narrow metallic contacts, depicted by the green 
atomic sites arranged in a square lattice, forming a normal- 
conductor/graphene/normal-conductor junction. The figure 
is not to scale, simulated devices were considerably bigger 
with ^ 700 carbon atoms along the zigzag edge. 
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is the annihilation (creation) operator acting on site i, 
the second sum is over nearest neighbors (i, j), and jij 
denotes the hopping matrix elements. The on-site poten- 
tial changes due to contributions from impurities and the 
gate voltage applied to the device, = e 



imp 



^gate- 



The 



graphene section consists of atomic sites placed in the 
hexagonal lattice with lattice constant a. As outlined 
in Ref. [32 , metallic contacts are modeled as a region 
with a square lattice arrangement with lattice constant 
dL = V^a, matching the A(B)-atom in the zigzag inter- 
face at the left (right) of the graphene flake. 

Disorder is introduced in two different ways. To gen- 
erate edge-disordered samples, we first consider surface 
disorder samples produced by randomly removing 30% of 
sites from the three outer atomic layers [20] . The second 
mechanism generates bulk disorder in the graphene de- 
first by choosing Nimp lattice sites at random 



vice I 

from the total number of sites in the device, with the 
ratio kept as N imp /N tot = 0.022. For the chosen lattice 
sites located at i?^, the on-site potential amplitude Vi is 
chosen randomly in the interval {—5^5) and smoothed out 
over a range ^ by convolution with a Gaussian, 



Ni 



^zmp 



mp 



i=l 



Vicxp 



Ri\ 



2^2 



(2) 



The range of the convolution is important given that in 
the extreme case of ^ << a the atomic scale disorder 
potential would break the A-B symmetry caused by hav- 
ing two atoms in the unit cell. For our calculations we 
assume that ^ = 2a resulting in a short-range potential 
(when compared to other length scales in the system) 
that varies smoothly on the atomic scale and suppress- 
ing the effect of intervalley scattering. 

The NEGF formalism [STJ |34] is a sophisticated frame- 
work for obtaining the transmission and other quantities 
in realistic devices. The retarded Green function in the 
atomic orbital basis is given by 

G^{E) = [E-H-^{E)]-\ (3) 

where the non-Hermitian self-energy matrix 

^{E)=^UE)+^r{E), (4) 

introduces the effect of attaching leads to the left and 
right ends of the device. The energy E tunes the Fermi 
level from the Dirac point {E = 0) to any charge den- 
sity induced by the gate voltage in the device. The self- 
energies determine the escape rates of electrons from the 
device into the semi-infinite ideal leads. Using the Lan- 
dauer formula ^35^, it is possible to obtain the transmis- 
sion function 



where G^{E) = {G^{E))'^ is the advanced Green func- 
tion. The matrices 



Tp{E, Vds 



E 



eVds 



E 



eVds 



(6) 

introduce a level broadening due to the coupling of the 
leads and a source-to-drain voltage given by Vds- Finally, 
the current in 2-terminal devices can be obtained from 
the Landauer formula 



I{Vsd) 



2e 



/oo 
dET{E) [f{E - iil) - f{E - m)] 
-CO 



T{E, Vds) = Tr [Tr{E, Vds)G^TL{E, 



(5) 



(7) 

where the energy window is defined from the difference of 
the Fermi functions of the macroscopic reservoirs where 
the leads terminate. In our calculations we assume the 
linear response limit (V^^ 0), where the relationship 
between conductance G and current is given by / = GVds • 



IV. RESULTS AND DISCUSSION 
A. Transport and the probability density function 

We investigate rectangular graphene sheets connected 
to ballistic, square lattice source and drain leads. As il- 
lustrated in Fig. [l] for all the results presented here, the 
edges of the graphene along the transport direction are in 
the armchair configuration. In order to extract the PDF 
for the charge flow, we calculate the conductance as a 
function of the position of the drain lead with the source 
lead fixed. The resulting quantity G{Ay), where Ay is 
the displacement of the drain lead, can be seen as an ap- 
proximation of the PDF (while it contains higher order 
effects such as back-reflections from the drain lead which 
is not completely transparent). We refer to G{Ay) as 
the displacement-conductance. The displacement of one 
of the leads from the center of the device is a theoreti- 
cal construct that can be easily investigated numerically 
but where the experimental implementation is not trivial. 
Yet, 'moving the leads' could be realized by fabricating 
several samples where the lead is a attached at a differ- 
ent positions for each device. Another possibility for the 
experimental realization is to replace one of the metallic 
leads by a scanning tunneling microscope tip. When the 
tip is brought close enough to the edge, the tunneling 
current could be measured as a function of tip position. 

In Fig. |2j we present the conductivity as a function of 
the displacement of the drain lead Ay along the edge of 
the sample for several values of the length of the system 
(L=ll, 21, 43 nm) and disorder models. When the two 
leads are collinear, the value of the displacement is {Ay = 
0). In the calculations presented in Fig.[2j both the width 
of the device and the width of the leads are kept constant 
at W = 74 nm and W^ = 10 nm, respectively. Thus, 
the schematic shown in Fig. [l] is not representative of 
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FIG. 2: Conductance versus displacement of one of the leads 
from the center of the device for Fermi energy Ef = O.57 
(with hopping matrix element 7). Each of the curves corre- 
sponds to different length of the system (L = 11, 21, and 43 
nm) while the width of the graphene device and of the leads 
are kept constant at W = 74 nm and = 10 nm, respec- 
tively, (a) Clean system, (b) System with edge disorder, (c) 
System with bulk disorder. 

the size of the devices considered, which are considerably 
bigger with about 700 atomic sites along the zigzag edge. 
To handle these system sizes we make use of efficient 
recursive Green's function methods. 

We explore how the conductance profiles behave as a 
function of the disorder model present by exploring three 
different systems. Fig. |2] (a) corresponds to a pristine 
graphene device, while (b) corresponds to a realization 
with edge disorder where 30% of the atoms at the outer 
edges are missing, and (c) to a system in the presence 
of disorder that would form electron-rich and hole-rich 
puddles [36 . Fig. [2] (a) shows that the maximum con- 
ductance is associated with the two leads being collinear 
{Ay ^ 0). For the clean device the maximum for the 
shorter device is at point Ay = 0, while for the other 
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FIG. 3: The transport mechanism is classified according to 
the shape of the central peak (of width 2Wl centered around 
Ay=0) of the conductance curves G{Ay) for (a,b) pristine 
system (c) system with edge disorder and (d) system with 
bulk disorder. (a,b) The central peak exhibits a triangular 
envelope indicating ballistic transport (while the Gaussian fit 
in (b) for the same conductance curve as in (a) does not de- 
scribe the cusp at A^=0 accurately). (c,d) The envelopes of 
the central peaks are well described by a Gaussian indicat- 
ing diffusive transport. Other parameters are Ep — 0.4997, 
L = 11 nm, W = 74 nm and Wl = 10 nm. 

two lengths there are local maxima near this point. Yet, 
there is an envelope of maximum conductance that is 
found when the vertical positions of the two leads over- 
lap with each other. Increasing the length of the device 
reduces the conductance peak due to the increased num- 
ber of reflection at the boundary of the sample. Similar 
features are seen in the curves for disordered samples. We 
find that for all three models (no disorder, edge disorder, 
bulk disorder) the width of the central peak is 2Wl' 

The characteristics of the transport can be studied via 
the conductance profiles. The profiles of the transmit- 
ted intensity are plotted in Fig. [3] with fits to either a 
linear or Gaussian curves. Diffusive transport is based 
on the multiple scattering paths taken by charge carriers 
as they transverse a device. In the case of conventional 
ballistic transport, the classical expectation of the PDF 
of the charge flow at the edge is a triangular shape of 
width 2Wl' That means, the triangular shape is the re- 
sult of a convolution of two square windows, where the 
window sizes are given by the lead width Wl- This dis- 
tribution for a device of length L = 11 nm is shown in 
Fig. [sj^a) where the envelope is well described by a tri- 
angular curve. In Fig. [sj^b), we show that near the cusp 
of the curve (at displacement zero) a Gaussian envelope 
differs significantly from the conductance curve. We em- 
phasize that a similar approach to distinguish normal 
diffusion (reflected by a Gaussian profile) from ballistic 
and Levy transport was made in [30]. In contrast, the 
conductance profiles for disordered systems are well de- 
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scribed by a Gaussian, as seen in Fig. [sj^c) for an edge 
disordered system and in Fig. [sj^d) for a system with bulk 
disorder. 

While the envelope of the central peak of width 
2Wl can be described in a quasiclassical way, the 
displacement-conductance curves exhibit several quan- 
tum interference effects. The regularly spaced fluctua- 
tions seen in the profiles for the clean systems are due to 
Fabry-Perot conductance interferences. Also, in Fig. |2] 
we observe fluctuations in the tails of the conductance 
curves. We carefully checked that these fluctuations are 
not only present for individual conductance curves and 
remain after averaging over an energy window. Thus, the 
fluctuations are consistent with universal conductance 
fluctuations (UCF) [371 EH] ^ which we therefore review. 

In a disordered mesoscopic conductor where the sample 
is smaller than the phase coherence length of the charge 
carriers but large compared to the average impurity spac- 
ing, the transmission is the result of the interference of 
many different paths through the system. As these paths 
are typically long compared to the wave length of the 
charge carriers, the accumulated phase along the paths 
changes randomly when an external parameter is varied 
(e.g. the magnetic field or gate voltage). This results in 
a random interference pattern, i.e. reproducible fluctua- 
tions in the conductance of a universal magnitude on the 
order of 2e^/h [371 [38]. Universal conductance fluctua- 
tions can also be created by the displacement of a single 
scatterer [391441] . But how can UCF occur in a ballistic 
device? The role of disorder in providing a distribution of 
random phases can as well be taken by chaos. Thus, bal- 
listic mesoscopic cavities like quantum dots in high mo- 
bility two-dimensional electron gases that form chaotic 
billiards show the same universal fluctuations (42fl44j . 
Strikingly, experiments in graphene quantum dots found 
strong indications for a chaotic billiard [45] . 

B. Charge flow 

Scanning probe microscopy (SPM) [16] has recently 
been used to image mesoscopic transport effects such 
universal conductance fluctuations [46], weak localiza- 
tion [47], and observed localized states [48| in graphene 
devices. The tip of the cryogenic SPM capacitively cou- 
ples to the graphene device inducing a movable scatterer. 
Conductance maps created via SPM have been proposed 
as a way to probe the chiral nature of the charge carriers 
in graphene . The work carried out by Berezovsky et 
al. [46l |47] and Schnez et al. [48 probed the regime of 
coherent diffusive transport in graphene. 

In the following, we discuss the possibility of using 
SPM to gain insight into the transport mechanisms and 
ballistic to diffusive cross-over in graphene. We simulated 
the effect of the capacitively coupled tip of the SPM to 
the graphene device as a point charge q located a height 



a above the substrate inducing a local perturbation to 
the charge density 

—qa . . 

where p is the in-plane radial coordinate away from the 
position of the point charge. The charge q is chosen to 
yield an rms charge density n ~ 4 x 10^^ e cm~^, in agree- 
ment with the observed charge puddles in graphene. The 
tip height a controls the width of the induced charge den- 
sity; we have chosen a = 10 nm, a distance from the 
substrate that is routinely used in SPM experiments and 
resulting in a half-width of the induced density of about 
20 nm. A conductance image is simulated by adding to 
the on-site potential the perturbation created by the tip 
to the contributions as well as the effect from the gate 
voltage and impurities present, = Cimp + ^gate + ^Up- 

In Fig.[4|we plot the conductance due to the presence of 
the tip as a function of the tip position. In both images, 
the effect of narrow lead constrictions are seen from the 
band of near constant conductance seen near the middle 
of the device. The simulated images show conductance 
fluctuations of order ~ ±2e^//i. The analytical results 
of Braun et al. [49 consider ballistic trajectories from 
the source to the drain leads in the presence of the tip 
potential. To first order, as treated in Ref. [49], there 
are only two trajectories that interfere with each other, 
source-to-drain trajectory and the source-tip potential- 
drain trajectory, revealing a spatial pattern due to the 
interference. Our numerical calculations are able to test 
the analytical results and give us insight into the robust- 
ness of the interference pattern to disorder while treating 
all possible trajectories from source to drain leads. 

The results presented in Fig. [4] calls for two important 
comments. First, for pristine devices, the interference 
pattern is present, with a clear suppression of the conduc- 
tance along bands parallel to the transport. The effect of 
multiple scattering reduces the interference pattern near 
edge of the device. Second, our calculations show an im- 
portant difference for disordered systems (right panel in 
Fig. [4]). In this case, for a device with high mobility and 
which would exhibit a long mean free path, the interfer- 
ence pattern is lost. In the case of pristine graphene the 
mean free path is no longer defined. However, even a 
small amount of disorder would make for some distribu- 
tion of length steps so that the mean free path can be 
defined, even if it still of a size comparable to the device. 

Finally, we would like to comment on possible experi- 
mental realizations. The SPM measurements have been 
carried out in current experimental setups for graphene 
devices on a Si02 substrate [46-48 and can be readily 
carried out in suspended graphene, where the ballistic 
transport regime can be reached, and in devices with 
hexagonal boron nitrate as substrate, where the ballistic 
to diffusive cross-over should be observable. 
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FIG. 4: SPM conductance maps. G{x,y) as a function of tip position, (left) pristine graphene sample (right) disordered 
device. The mobility of the disordered device is comparable to that of devices in crystalline substrates and the mean free path 
comparable to system size. 



V. CONCLUSION 

We proposed a method to study the mechanisms of 
electronic transport in graphene that is based upon dis- 
placement of the drain lead along the zigzag edge. The 
method reflects the PDF of the charge flow at the edge of 
the device. Several features based upon quantum inter- 
ference effects (Fabry-Perot resonances and UCF) were 
found in the resulting conductance curves. The gen- 
eral method that we applied to graphene can be used as 
well to study other systems such as semiconductor het- 
erostructures. Moreover the method renders a distinc- 
tion between normal and anomalous transport possible. 
In the second part of our study, we simulated SPM mea- 
surements for the same devices (but with flxed leads), 
which illuminates the flow of charge inside the device 
and complements the transport calculations. 
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